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Interference experiments with independent condensates provide a powerful tool for analyzing 
correlation functions. Scaling of the average fringe contrast with the system size is determined by 
the two-point correlation function and can be used to study the Luttinger liquid liquid behavior 
in one-dimensional systems and to observe the Kosterlitz-Thouless transition in two-dimensional 
quasicondensates. Additionally, higher moments of the fringe contrast can be used to determine 
the higher order correlation functions. In this article we focus on interference experiments with 
one-dimensional Bose liquids and show that methods of conformal field theory can be applied to 
calculate the full quantum distribution function of the fringe contrast. 

Correlation functions provide a convenient language for characterizing quantum states. First-order coherence in 
optics, which underlies classical Young's interference experiments, is mathematically equivalent to the factorizability 
of the first order correlation function (xi, x% ) = (xi )E(+) (x 2 )), where .EW are the electric field components 

which vary as e ±luJt . Photon bunching in Hanbury-Brown and Twiss experiments can be understood using the second 
order correlation function G^{t) = (£?(") {t)E { -^ [t + t)E { +^ (t + t)E {+ ^ (t)> . The same correlation function describes 
photon antibunching in quantum mechanical states of light such as number states of photons [lj]. Analogously, in 
condensed matter physics many common experiments can be understood as probes of the appropriate two point 
correlation functions. For example, conductivity measurements (dc or finite frequency) correspond to current-current 
correlation functions [2], angle resolved photoemission 3, 4] and tunneling experiments 0] probe single particle Green's 
functions, X-Ray and neutron scattering experiments can be related to charge and spin correlation functions Q. 

An important question for current experiments with ultracold atoms is finding new methods of characterizing 
strongly correlated many-body states, such as systems near Feshbach resonances, rotating condensates, atoms in 
optical lattices, low-dimensional systems (see ref. Q for a review). The existing experimental toolbox includes 
Bragg scattering, which measures the dynamic structure factor [8f (i.e. the imaginary part of the density-density 
correlation function), and the RF spectroscopy (see e.g. Refs. 0, [l(| ), which is essentially equivalent to finite 
frequency conductivity measurements in solid state systems. These techniques measure different types of two point 
correlation functions. An interesting question to ask is whether one can do experiments with systems of cold atoms 
that would measure higher order correlation functions. In two recent papers fill [l2j we argued that this can be done 
using interference experiments with independent condensates (see also Ref. 13]). In this article we will review these 
ideas while focusing on a specific case of one-dimensional condensates. 

We mention in passing that another interesting technique, which is unique to systems of cold atoms, is the time of 
flight experiments (see e.g. Ref. [HI])- When expansion of atoms released from the trap is ballistic, the momentum 
distribution of atoms inside the original system gets mapped into the density distribution after the expansion. Mea- 
surements of quantum noise after the expansion have been used to study Hanbury-Brown- Twiss correlations for bosons 
and fermions in optical Iattices[l5|,[l(|[l7j and pairing correlations on the molecular side of the Feshbach resonance [l8j . 
The correlation function measured in such experiments is a non-local operator from the point of view of the original 
system but should provide "smoking gun" signatures of several many-body phases including antiferromagnetically 
ordered Mott states and paired states of fermions [l9l|. 

The interference experiments we consider are shown in Fig. [TJ Two independent quasi-condensates are allowed to 
expand in the transverse direction. After sufficient expansion, the atom distribution is measured using the imaging 
beam. Everywhere in this paper we consider the two condensates to be identical, although our analysis can be 
generalized to the case of different condensates. In the set-up we discuss, the measuring beam is oriented along the 
axis of the original cloud. This beam integrates over local interference patterns and the fringes that we see on a 
detector are a result of averaging over the entire imaging length. In the presence of thermal or quantum fluctuations, 
interference patterns are not in phase at different points and the resulting interference fringes have a reduced contrast. 
Such reduction carries information about fluctuations in the original clouds. The main topic of this article is a non- 
trivial information that we can extract by analyzing such reduced contrast. We point out that an alternative set-up, 
which is also possible with current experiments, is to measure not along the system axis but from the side. This 
provides a local picture of interference fringes rather than the integrated one. A convenient way to analyze such 
data would be to integrate interference patterns numerically and study how the average contrast changes with the 
system size (see discussion below) . So the main advantage of the latter set-up is a dramatic reduction in the number 
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FIG. 1: Experimental setup 

of measurements needed to determine the scaling of interference fringes. Conceptually, however, analysis is identical 
to the case that we discuss in this paper. We assume that before the expansion, atoms are confined to the lowest 
transverse channels of their respective traps and that the optical imaging length L (which is smaller than the size of 
the system in the axial direction) is much larger than the coherence length of the condensates. This allows us to use 
an effective Luttinger liquid description of the interacting bosons (20|. 

The quantum observable corresponding to the interference amplitude of the two condensates is given by Aq — 

Jq dza\(z)a2(z). Here a\ and a 2 are the bosonic operators in the two systems before the expansion, and the integrals 
are taken along the condensates. L is the imaging length that is in general smaller than the full condensate length. If 
the two condensates are decoupled from each other, the expectation value of {Aq) vanishes. This does not mean that 
\Aq\ is zero in each individual measurement but rather shows that the phase of Aq is random. Hence the position 
of interference fringes is completely random from shot to shot. The quantum mechanical expectation value is defined 
as a result of averaging over many experimental runs. Superimposed interference patterns with random phases wash 
each other out and interference fringes averaged over many shots disappear [IJ. However, what we are interested in 
is the amplitude of the fringes in an individual measurement. Thus we should consider an observable that does not 
involve the random phase of Aq. We take 

\ A q\ 2 = / dz 1 dz2a\(zi)a 1 (z2)a\(z2)a 2 {zi) (1) 



Jo Jo 

To simplify calculations we consider the limit of a large system when we are allowed to take the normal ordered 
expression 22]. Taking the expectation value of the last equation produces a nontrivial result since on the right hand 
side operators that correspond to different clouds decouple. Assuming that the condensates are identical gives 

{\Aq\ 2 )=L [ L dz(aHzHO)) 2 - (2) 
Jo 

Thus the average intensity of the interference pattern depends on the two point correlation function along the indi- 
vidual one-dimensional condensates. In this regime the long distance correlations decay as a power law within the 
imaging length L, (at(z)a(O)) ~ z~ 1 ^ 2K , where K is the Luttinger parameter. For bosons with a repulsive short-range 
potential, K ranges between 1 and oo, with K = 1 corresponding to strong interactions, or "impenetrable" bosons, 
while K — > oo for non- interacting bosons. Using the power law correlations in ([2]), we find that the interference inten- 
sity scales as a non-trivial power of the imaging length (|Aq| 2 ) = L 2 ~ x / K . This property can be used as a sensitive 
probe of the Luttinger parameter. In the non-interacting limit (K — > oo), the scaling is linear, Aq = y/ (\Aq | 2 ) ~ L, as 
expected for a fully coherent system. Interestingly, in the impenetrable limit, Aq ~ y/~L, as for short range exponen- 
tially decaying correlations. One may ask whether it makes sense to discuss zero temperature limit when experiments 
are always done at finite temperature. It is important to realize that experiments are done in systems of finite size. 
At finite temperatures and at long distances, correlations always decay exponentially. However, long distances means 
distances larger than the correlation length £y ~ v s /T, where v s is the sound velocity. As the temperature is lowered, 
£t increases and at some point becomes larger than the system size. At this point (and at all lower temperatures), 
correlations within the system size are those of the zero temperature system. It may be useful to make one remark 
regarding experimental constrains for observing quantum Luttinger liquid behavior with ultracold atoms. A common 
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parameter used to characterize the strength of interactions between particles is the ratio of the interaction energy 
to the kinetic energy 7 = mg/Ti 2 n. Here m is the mass of atoms, g is the interaction strength, and n is the one- 
dimensional density. To observe strongly interacting Bose liquid we want 7 to be large, which may be achieved by 
reducing the density. However we also need to satisfy the condition £t 3> L. If we use Bogoliubov expression for v s 
we find L Wng/m >> T. Hence for a smaller density it takes a lower temperature to reach the quantum limit. In 
experiments one needs to find the balance between having the interesting regime of strong interactions and reaching 
the quantum limit. 

We point out that analysis presented above for one dimensional systems can be easily generalized to interference 



between planar condensates at finite temperatures 11]. The power law decay of correlations, and thus the scaling of 
the interference intensity with imaging size is related in this case to the supcrfluid stiffness of the condensates. This 
was used by Hadzibabic et al [23| to observe the universal jump of superfluid stiffness across the Kosterlitz-Thouless 
transition in pancake condensates. 

So far we have discussed the average intensity of the interference pattern and showed that it contains information 
about two-point correlation functions. It is important to realize that interference experiments correspond to a classical 
measurement of a quantum mechanical state. Hence they contain an intrinsic quantum mechanical noise and the result 
of each individual measurement will be different from the average value. Generalization of the argument that led to 
eq. ([2]) shows that higher moments of the distribution function of interference amplitudes correspond to high-order 
correlation functions. Hence the knowledge of the entire distribution function reveals global properties of the system 
that depend on very high order non-local correlation functions. 

From equation |T]) one finds that for two one-dimensional condensates higher moments of the interference fringe 
amplitude are given by (Tl| 



(\A Q n = AfZ 2n , where A = y/ Cp& /K L^l* , (3) 

where C is a constant of order unity, p is the particle density in each condensate, £/, is the short range cutoff equal to 
the healing length. The coefficients Zi n in Eq. Q are given by[24j: 



z 2n (K)= / .../ n 
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The coefficients Zi n originally appeared in the grand canonical partition function of a neutral two-component 
Coulomb gas on a circle 

Z ( K ^) = Ht^ Z ^ K )- ( 5 ) 

n=0 ^ '' 

Here x is the fugacity of Coulomb charges and Zi n describes contributions from configurations with 2n charges (i.e. 
canonical partition functions). The partition function ([5|) with K > 1 describes several problems in statistical physics 
(see Ref [25| and references therein). It is also related to the problems of an impurity in a Luttinger liquid fill. |26| 
and Kondo effect^. 

When describing interference experiments it is convenient to define the normalized amplitude of interference fringes 
a = \Aq\ 2 /Aq. From ^ we find that (a 2n ) = Z 2n , so by performing experiments that measure the distribution 
function W(a), we get direct access to the partition function ([5]). We point out that W(a) can be used to compute all 
moments of |^4q| 2 , and therefore contains information about high order correlation functions of the interacting Bose 
liquids. 

From the Taylor expansion of the modified Bessel function we find 

/>oo 

Z(K,x)= W{a)I {2xy/^)da. (6) 



Inverting Eq. © allows one to express the probability VT(a) through the partition function Z(K,x). Noting that 
Io(ix) = Jo(x) and using the completeness relation for Bessel functions, J °° Jo(Xx) Jo(Xy)\x\\d\ = S(\x\ — |y|), we 
obtain 

poo 

W(a) = 2 Z(K,ix)J (2xVa)xdx. (7) 
Jo 

The remaining problem is to calculate the partition function Z(K, ix). In this paper we will use a method based 
on the studies of the integrable structure of conformal field theories [28] . In particular, it was shown that the vacuum 
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expectation value of Baxter's Q operator (central to the integrable structure of the models [29|]), coincides with the 
grand partition function of interest: 



3 (X) = Z(K,-ix), 



(8) 



where x is related to the spectral parameter X, x = nX/ sm(n/2K). It was conjectured in Refs. [3Cj,|31| that the vacuum 
expectation value Q vac (A) is proportional to the spectral determinant of the single particle Schrodinger equation 
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(9) 



So 



g vac (A) = D(pX 2 ), where p = {AKf- x l K [T(l - 1/(2K))} 2 , D(E) is the spectral determinant defined as D{E) 
x (l — E/E n ), and E n are the eigenvalues of ([9]). Thus, we have 



Z(K, ix) = JJ ( 1 

n=l 



E n 



(10) 



The distribution function W(a) is shown in Fig. [2]for several values oiK. For K close to 1 (Tonks-Girardeau limit), 
W(a) is a wide Poissonian function, which gradually narrows as K increases, finally becoming a narrow 5-function at 
K — > oo (limit of noninteracting bosons). Note that the distribution function remains asymmetric for arbitrarily large 
K. In fact we find that W(a — 1), where a — a/ (a) = \A\ 2 /(\A\ 2 }, tends to a universal scaling form, parameterized 
by a single number characterizing the width of the distribution: 5 = (a 2 ) — 1 it j\f%K. We conjecture that the 
limiting form of is the Gumbel distribution function [l2j]. The appearance of the Gumbel distribution in this 

problem is not surprising. This distribution was introduced to describe rare events, such as earthquakes or stock 
market crashes, which act in one direction. For example, earthquakes destroy property but do not create it. In the 
limit of large K, bosonic quasicondensates have small fluctuations and exhibit good interference patterns in most 
cases. But occasionally there are strong fluctuations which lead to an appreciable decrease of the contrast. This gives 
rise to a strong asymmetry of the resulting distribution function. 
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FIG. 2: Evolution of the distribution function W(a) for different values of K. 
the delta- function (see the text). 



At larger values of K the function W tends to 



The construction of [28( is based on the representation space of the Virasoro algebra for the set of central charges 
satisfying c = 1 - 6{V2K - 1/V2K)) 2 and the highest weight A = (c - l)/24 which gives c < -2 for K > 1. For all 
K > 1 the central charge of the Virasoro algebra is negative. Theories with negative central charges appear in different 
contexts in statistical mechanics, stochastic growth models, 2D quantum gravity, models of 2D turbulence and even 
high-energy QCD. The experimentally measured distribution function W can be inverted (using the Bessel functions 
completeness relation) to obtain the Q-operator which in turn can be used to reconstruct the transfer matrices of 
the above mentioned models with negative c. Therefore, the interference of condensates provides a possible way to 
explore the interesting physics of various models ranging from statistical to high energy physics. 

To summarize, in this paper we discussed interference experiments between two independent one-dimensional quasi- 
condensates. We showed that quantum phase fluctuations act to suppress the average interference contrast and to 
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FIG. 3: Limit of large K. Scaled distribution function SaW ((a — l)/8a), where a = a/ '{a) = |j4| 2 /(|A| 2 ) and 5a is the width 
of the distribution. The function W is multiplied by 5a to preserve the total probability, which must be equal to unity. The 
dashed and dotted lines correspond to different values of K . The solid line corresponds to the conjectured Gumbel distribution. 



induce fluctuations in this quantity from one experimental run to the next. The average interference contrast depends 
on the two point correlation. It therefore scales as a non trivial power of the imaging length related to the power 
law decay of the correlations. The distribution function characterizing the shot to shot fluctuations of the contrast 
contains information on high order correlation functions. We computed the distribution function of the amplitude of 
interference fringes using the relation of this problem to the partition function of the logarithmic Coulomb gas on a 
circle and to the properties of Q operators of conformal field theories with negative central charges. We showed that 
the distribution function of fringe amplitudes is related to non-trivial physical properties of a variety of interesting 
statistical and field-theoretical models. 
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